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Abstract 

A /i-uniform hypergraph H — (V, E) is called (£, fc)-orientable if there exists an assignment of 
each hyperedge e £ E to exactly £ of its vertices v G e such that no vertex is assigned more than k 
hyperedges. Let H njm , t h be a hypergraph, drawn uniformly at random from the set of all ^.-uniform 
hypergraphs with n vertices and m edges. In this paper, we determine the threshold of the existence 
of a (I, fc)-orientation of -ffn.m.h for k > 1 and h > £ > 1, extending recent results motivated by 
applications such as cuckoo hashing or load balancing with guaranteed maximum load. Our proof 
combines the local weak convergence of sparse graphs and a careful analysis of a Gibbs measure on 
spanning subgraphs with degree constraints. It allows us to deal with a much broader class than the 
uniform hypergraphs. 
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1 Introduction 

Motivated by load balancing problems [37], Gao and Wormald [IT introduced the following generalisation 
to random hypergraphs of a commonly studied orientation problem on graphs. A /i-uniform hypergraph 
is a hypergraph such that all its hyperedges have size h. Let h > £ be two given positive integers. A 
hyperedge is said to be ^-oriented if exactly £ distinct vertices in it are marked with positive signs with 
respect to the hyperedge. The indegree of a vertex is the number of positive signs it receives. Let k be a 
positive integer. A (£, k) -orientation of an /i-uniform hypergraph is a ^-orientation of all hyperedges such 
that each vertex has indegree at most k. If such a (£, fc)-orientation exists, we say that the hypergraph 
is (£, fc)-orientable. We consider Gn,m,h the probability space of the set of all /i-uniform hypergraphs on 
n vertices and m hyperedges with the uniform distribution. A random /i-uniform hypergraph is then 
denoted by H n ^ m ,h- We are now ready to state our main result in this framework: 

Theorem 1. Let Q(x,y) — e~ x X^>y ff an< ^ Bin(n,p) denote a binomial random variable with param- 
eters n£N and p £ [0,1], i.e. F(Bin(n,p) = k) = (l)p k (l - p) n ~ k . For integers h> £>l, k>l with 
max(ft — £,k) > 2, let £* be the unique positive solution to 

hk = c E [max (£ - Bin(h, 1 - Q(C, k)),0)] 



Q(e, k + 1)P {Bm(h -1,1- Q(f, A)) < ' 
Let 

r 



c h,l,k 



h¥ (Bin(h -1,1- Q{£*,k)) < £) ' 
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Then 



lim P (H n i cn i h is (£, k)-orientable) = 

n— too ' 

The characterisation of the threshold c£ fc £ in [T7] (for fc sufficiently large) involves the solution of a 
differential equation system which is rather complicated (according to the authors themselves) and does 
not allow to get explicit values for c* h k t We believe that our method of proof and the characterisation 
of the threshold in Theorem [1] is much simpler. Note that for the case k = 1 and I = ft- — 1, the threshold 
for orientablity is equal to the threshold for the apparition of a giant component so that we have 

1 

Cfc . fc - 1 - 1 ft(ft-l)' 

Hence our result allows to compute the orientation threshold on the whole range of parameters k > 1 and 

£c* 

h > £ > 1. To illustrate our result, we computed numerical values of the critical load h k ' k for different 
values of the parameters. First when £ and k vary while ft = 4 is fixed: 



£/k 
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1 


0.9767701648 


0.9982414840 


0.9997951433 


0.9999720662 


0.9999958680 


0.9999993570 


2 


0.7596968140 


0.9266442602 


0.9676950000 


0.9834603210 


0.9908051880 


0.9946173050 


3 


0.0833333333 


0.6612827547 


0.7892143791 


0.8525202000 


0.8898186996 


0.9141344769 


Then when ft and k vary while £ = 2 is 


fixed: 








h/k 
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0.8833250296 


0.9730747564 


0.9909792334 


0.9964896324 


0.9985201920 


0.9993444714 


G 


0.9378552354 


0.9894605852 


0.9974188480 


0.9992698236 


0.9997769140 


0.9999284650 


7 


0.9652101902 


0.9957801256 


0.9992689074 


0.9998543770 


0.9999687056 


0.9999929390 



We should also stress that our result (and our proof) unifies various results available for different 
ranges of the parameters. 

The case £ = 1 has attracted a lot of attention and can be described in different terminologies. In 
the balanced allocation paradigm [3], we have m balls and n bins. To each ball, two bins are assigned at 
random. Each ball is to be placed in one of the two bins assigned to it; the aim is to keep the maximal 
load small. Another formulation of the same problem can be given in data structure language: in the 
cuckoo hashing method fM\ , each one of m keys is assigned two locations in a hash table of size n and 
can be stored in one of the two locations. If each location has capacity one (or in the balanced allocation, 
maximal load cannot exceed one) , the offline version of this problem corresponds exactly to the standard 
1-orientability of the classical random graph G n . m drawn uniformly from the set of all graphs with n 
vertices and m edges. To see the connection, associate to each location a vertex of the graph and to each 
key an edge of the graph: the orientation of the edges correspond to the allocation. It is easily seen in 
this case that we must have m < n/2 in order for G n ^ m to be 1-orientable. An interesting generalisation 
of this problem considers bins/locations of capacity k > 1 [IS]. This generalisation (with the number 
of choices per ball/key still equals to 2) corresponds to the (1, fc)-orientability described above for the 
random graph G n>m = H n ,m,i- For k > 2, the sharp threshold for the fc-orientability of the random 
graph Gn iTn , corresponding to c* 2 1 k in our Theorem [T] was found simultaneously by Cain, Sanders and 
Wormald [5] and Fernholz and Ramachandran [13]. Another generalisation allows for ft > 2 choices 
of bins/locations. In this case, the graph associated to the problem is a random ft-hypergraph. The 
(1, l)-orientability threshold of ft-uniform random hypergraph with ft > 3, corresponding to our c* h 1 1 in 
Theorem [1] has been independently computed by Dietzfelbinger et al. [TT], Fountoulakis and Panagiotou 
[15] and Frieze and Melsted [16]. This corresponds to a case where bins/locations have unit capacity. 
Recently, the extension to capacity k > 1 has been solved by Fountoulakis, Kosha and Panagiotou [Tl] 
for any ft > 3 (refining results of [T7] for this particular case). Our derivation of c* h 1 k agrees with [T3]. In 
the generalisation proposed in [T7] , each batch of £ balls has ft > I choices of bins and each of the n bins 
has a capacity k. Our Theorem [T] gives the threshold Ch,i,k such that, as n goes to infinity: if n/m < c 



| ifc> c* hek , 
\ 1 ifc<c*. 
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with c < Ch.e.k then there exists with high probability an allocation of the m balls such that the maximal 
capacity of a bin is one; and if n/m > c with c > Ch.e.k then such allocation does not exist with high 
probability. 

Our approach is completely different form the works cited above. We consider the incidence graph of 
the /i-uniform hypergraph H, which is a bipartite graph G = (AUB,E) where A is the set of hyperedges, 
i.e. are vertices in G with degree h and B is the set of vertices of H . We then consider spanning subgraphs 
S = (A U B, F) of G with degree constraints: any vertex from A must have degree at most £ in S while 
any vertex from B must have degree at most k in S. Let the size of such a spanning subgraph be the 
number of edges |F| in S. The following claim is easy to check and will be the basis of our approach: 
H is (£, fc)-orientable if and only if all vertices in A have degree £ in any maximum spanning subgraph 
with degree constraints (£,k). In this case, the size of any maximum spanning subgraph is £\A\. Indeed 
in the case (£, k) = (1, 1) a spanning subgraph with degree constraints (1, 1) is simply a matching of G. 
Based on this observation, Bordenave, Salez and the author already derived the value of c* h 1 1 in [? . The 
analysis of maximum spanning subgraphs with general degree constraints £ and k requires a significant 
extension of the results in [7]. Wagner [25] and, more closely related to our work, Salez [25] study the 
generating polynomial for spanning subgraphs with degree constraints. It follows from 28 that for a 
sequence of graphs whose size goes to infinity and having a random weak limit (see definition in the sequel 
or [5], [3]), the rescaled size of a maximum spanning subgraph converges. In the case, where the random 
weak limit is a Galton Watson tree, this limit is characterised in [26, and computed in the particular 
case of constant degree constraint. In this work, we are able to simplify the characterisation of [55] in 
our Proposition [5] and to connect it to a simple message-passing algorithm. It allows us to bypass the 
resolution of a difficult recursive distributional equation (which was a key step in [7] or [26] ). It should 
perhaps be noted that the result in this paper is stronger than [5] [T3] [TTJ [T5J Q31 E] in the sense that it 
gives the size of the largest spanning subgraph for all values of the parameter c > 0. We state this result 
explicitly in the following theorem: 

Theorem 2. Denote by Mi^(H) the size of a maximum spanning subgraph of the bipartite graph H with 
degree constraints (£,k). With the same notation as in Theorem]]] and for any integers h > £ > 1 and 
k > 1, we define the function of (q,c) £ [0, 1] x R + : 

Tl, k {q, c) = E [min (£, Bin(h, 1 - Q{chq, k)))\ + k Q( ch ^ k + l ) 



Then we have for any c > 0, 



lim —M ek (H n i cn \ h )= inf F t , k {q,c) 

n-s-oc cn gG[0,l] 



It follows from calculations made in Section |5J that for c < c* l£fc , we have inf gg r 01 ] Ti tk (q, c) = £ 
while for c > e k , we have inf ge [ 01 ] J-"(. t k(q, c) < £■ In particular, for c > c^ fc , Theorem [5] gives 
the asymptotic fraction for the number of balls/keys that cannot be stored in the system. A similar 
statement was proved in [16] in the particular case of £ = k = 1 corresponding to spanning subgraphs 
being matchings. In the sequel, we will give a more general statement of this result (see Theorem[3]) which 
allows to deal with a larger family of random hypergraphs. Indeed, we believe that this approach will 
allow to deal with more complex situations where degree constraints could be random and/or asymptotic 
degree distributions could be changed and can lead to efficient load balancing algorithms. We refer to 
[21) where such extensions are explored. 

The rest of this paper is organised as follows. We give an overview of our proof in Section [5] We start 
with a careful analysis of a Boltzmann-Gibbs distribution on spanning subgraphs with degree constraints 
in Section [3] and show some crucial monotonicity of the model. We then give an explicit characterisation 
of the size of a maximum spanning subgraph when the underlying graph is a finite tree and are able to 
extend it to a possibly infinite tree using the important notion of unimodularity 2 . In Section [5] we 
apply these general results to the particular case where the underlying tree is a branching process. This 
allows us to derive the asymptotic for the size of a maximum spanning subgraph for a converging sequence 
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of bipartite graphs. In Section we apply our results for the particular sequence of graph H n .y C n\ji & n d 
derive Theorem [T] 



2 Gibbs measures and overview of the proof 

We consider a finite simple graph G = (V, E) with a vector of N v denoted by w = (w v , v £ V) and 
called the vector of (degree) constraints. We are interested in spanning subgraphs (V, F) with degree 
constraints given by the vector w. Each such subgraph is determined by its edge-set F C E encoded by 
the vector B = (B e , e £ E) e {0, 1} E defined by B e = 1 if and only if e £ F. We say that a spanning 
subgraph B satisfies the degree constraints or is admissible if for all v £ V, we have z^eedv^e — w v> 
where dv denotes the set of incident edges in G to v. We introduce the family of probability distributions 
on the set of admissible spanning subgraphs parametrised by a parameter z > 0: 

/4(B) = ^ (1) 

where Pg(z) — B " Huey l(S e edi> ^e — Wl >)- We a ^ so define the size of the spanning subgraph 

by \F\ = J2 e Be and denote the maximum size by M(G) = max{^ e B e : B admissible}. Those spanning 
subgraphs which achieve this maximum are called maximum spanning subgraphs. For any finite graph, 
when z tends to infinity, the distribution /j,q converges to the uniform distribution over maximum spanning 
subgraphs. For an admissible spanning subgraph, the degree of v in the subgraph is simply J2 e edv Be- 
By linearity of expectation, the mean degree of v under the law [x z G is D* := z~2 e ^dv (B e = 1) so that 
we have 

M(G) = ^ £ lim D *- ( 2 ) 

vev 



The main part of our work will be devoted to the computation of the limit on the right-hand side of 
@. In the remaining part of this section, we give an informal description of the main steps. The reader 
interested in the mathematical proof can skip the rest of this section and proceeds directly to Section [3J 

Recall from the Introduction that we see hypergraph H as bipartite graph G — (A U B,E). Then, 
an hypergraph is {I, fc)-orientable if and only if all vertices in A have degree £ in any maximum spanning 
subgraph of the corresponding bipartite graph with degree constraints (£, k). Indeed in this case, we 
have for any v € A, lim;,-^ D z v = £ so that the size of a maximum spanning subgraph is M(G) = £\A\. 
Our main result will show that liiUn^oo Ml^iMi — inf ? Jj^f^c), where G n (c) = (A n U B n ,E n ) is the 
bipartite graph associated to H n y cn ^ h and the function J^i^iQ^c) is defined in Theorem[5J In particular, 
when infq ^^(g, c) < £, our result allows us to conclude that the hypergraph is not (£, /c)-orientable for 
such values of c. In order to get the second half of the Theorem [TJ when va£ q J-^ t k{q, c) = £, we need to 
show that with high probability (as n tends to infinity) there are actually no vertices in A n with degree 
less than £ in a maximum spanning subgraph (since the limit ensures only that the number of such 
vertices is o(n)). For this part of the proof, we use a density argument which relies on a combinatorial 
argument of [T5]. This is done in Section [5J 

We now concentrate on the computation of the limit linin^oo using techniques from the 

objective method developed by Aldous and Steele [3]- A fundamental ingredient of the proof is the 
fact that the bipartite graphs associated to /i-uniform hypergraphs considered in this paper, i.e. with n 
vertices and cn hyperedges are locally tree-like: with high probability, there is no cycle in a ball (of fixed 
radius) around a vertex chosen at random. It is then instructive to study maximum spanning subgraphs 
when the underlying graph is a tree. Let first study the Gibbs measures defined by ([T]) in the limit z — > oo 
in order to analyse maximum spanning subgraphs. When the underlying graph is a finite tree, we can 
use a more direct and algorithmic way that we now describe. 
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To study the (t, fc)-orientability of the hypergraph H associated to the bipartite graph G, the vector 
of degree constraints w should be chosen such that w v — I for v E A and w v = k for v E B. For 
simplicity, we assume here that the vector of degree constraints is constant so that all vertices have the 
same degree constraint say w > 1. Consider now the following message-passing algorithm forwarding 
messages in {0, 1} on the oriented edges of the underlying tree G as follows: at each round, each oriented 
edge forwards a message, hence two messages are sent on each edge (one in each direction) at each round. 
The message passed on the oriented edge e = (u, v) is if the sum of the incoming messages to u from 
neighbours different from v in previous round is at least w and the message is 1 otherwise, i.e. if the sum 
of the incoming messages is strictly less than w. Let 1^ e {0, 1}-^ be the vector describing the messages 
sent on the oriented edges in E at the /c-th round of the algorithm. Denote by Vg the action of the 
algorithm on the messages in one round so that Ifc+i = Vc^-k)- Assume that the algorithm is initialised 
with all messages set to one: Iq = 1. Figure [T] shows an example with w = 2 for the messages exchanged 
for the three first rounds. In Figure [TJ iterating a fourth time the algorithm would again give I3. Indeed, 




it is easy to see that the algorithm will converge on any finite tree after a number of steps equals to 
at most the diameter of the tree, whatever the initial condition. Hence the messages of the algorithm 
converge to a vector I* = if G i^) solving the fixed-point equation I* = T>g(J*) and the size of a 
maximum spanning subgraph is given by 

^e(We j*>»+i)+i(e E^V ( 3 ) 

where dv is the set of oriented edges toward v. For example, one can check on Figure [1] using I3 = I* 
that the formula given by §3§ equals 5 which is the size of a maximum spanning subgraph. 

Proving the correctness of the algorithm and of the formula ([3]) on finite trees is simple once the 
following observations have been made: 

• a vertex v such that X^eSu ^i? — w ^> w ^ have degree w in any maximum spanning subgraph 
and all messages in I* sent by v will be zero. 

• an edge (u, v) such that = I£-y u — 1 will be covered (i.e. with -B( Uj „) = 1) by any maximum 
spanning subgraph. 
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• an edge (u, v) such that = Iv^. u = w ^ never be covered (i.e. with B( UjV \ ~ 0) in any 

maximum spanning subgraph. 

Note that the correctness of the algorithm is ensured for trees only, but the definition of the algorithm 
does not require the graph to be a tree. It makes only local computations and can be used on any graph 
(without guarantee of converging and even in this case without guarantee of correctness). Since the 
bipartite graphs associated to H n .\cn\,h are not trees but are locally tree like, it is tempting to use the 
algorithm directly on these graphs. It turns out that for low values of c, the algorithm will converge and 
will also be correct (with high probability). This is not a surprise since for c < 1/h, the random graph 
is essentially a collection of small tree components. It turns out that the algorithm allows to compute 
the size of a maximum spanning subgraph for values of c above 1/h but it breaks down at some higher 
value of c, indeed exactly when the (£, k + l)-core as defined in [T7] appears. We did not try to make this 
claim rigorous as it is not required for our analysis. In the particular case where h = 2 and I = k = 1 
(i.e. w v = 1 for all vertices v), our problem reduces to the problem of maximum matching in Erdos-Renyi 
random graphs and our claim follows from the analysis of Karp and Sipser [19] . More precisely, the greedy 
algorithm described above corresponds in this case to the standard leaf-removal algorithm studied in [115] 
(see also Proposition 12 in [5] and the Appendix of [7]). The (1, l)-core is simply called the core in these 
references and it appears if the mean degree is above e. If there exists a (£, k + l)-core, the intuition is as 
follows: if we approximate the graph by a branching process and run the algorithm on this tree starting 
from level k from a root, the influence of the boundary conditions on the value given by the algorithm 
at the root is positive as we let k tends to infinity. This translates into the fact that on an infinite tree, 
there might exist several solutions to the fixed-point equation I = "Pg(I)- It is then natural to ask which 
one is associated to the large n limit maximum spanning subgraphs. From an algorithmic viewpoint, 
there is 'no correlation decay' and the computations made by the algorithm is not anymore local. 

In order to bypass this absence of 'correlation decay', we borrow ideas from statistical physics by 
introducing the Gibbs measures /i G defined in ((T|) parametrised by a parameter z > (usually called the 
activity or the fugacity) |10j . Informally, the introduction of this parameter z will allow us to capture 
sufficient additional information on our problem in order to identify the 'right' solution to the fixed-point 
equation I = Vg(1)-, when we let z goes to infinity. Our first step in the analysis of these measures is to 
derive a message-passing algorithm allowing to compute the mean degree D% of any vertex v in a spanning 
subgraph taken at random according to the probability distribution /i G . We will proceed by first defining 
the local computations required at each node and we call them the local operators. We use these building 
blocks to define a message-passing algorithm which is valid on any finite tree. In particular, we show that 
as z tends to infinity, the dynamic of the algorithm becomes exactly the one described previously in this 
section. As a by-product, we prove the validity of ©, see Proposition [3] The advantage of considering 
these algorithms with z < oo is that we are able to define them properly on infinite graphs if these 
graphs have a natural stationarity property called unimodularity [5] . Section 13.31 presents this notion 
in details and shows how it is used in our framework. Proposition [4] shows that 'the message-passing 
algorithm converges to a unique fixed point for any z < oo'. In other words, Proposition U shows that 
for any z < oo, there is 'correlation decay'. In the limit z = oo, the algorithm may 'have more than 
one fixed-point' (corresponding to the absence of correlation decay mentioned above), but Proposition [5] 
allows to select the 'valid' fixed-point which gives the mean degree of a vertex 'picked at random' in a 
'maximum' spanning subgraph of a possibly infinite graph. The notion of unimodularity is a key concept 
in making these statements rigorous. This is the price to pay in order to work directly with the infinite, 
probabilistic object obtained as the limit of our finite problem as n tends to infinity. The reward of this 
objective method [3] is that in our case, this infinite object with fixed distributional properties is simple 
to analyse (see Section [5]) and captures all the necessary information on the asymptotic behaviour of the 
original sequence in order to compute the limit ([2]) . 
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3 Analysis of Gibbs measures 



This section is devoted to the computation of the limit appearing on the right-hand side of @. We first 
deal with the particular case where G is a star (i.e. a tree with one internal node and several leaves) 
then we show that the computation can be done recursively on a finite tree. We also show how these 
recursions can be extended to possibly infinite unimodular networks (defined in the sequel). 



3.1 Local operators 

In this section, we define local operators associated to the degree constraints defining an admissible 
spanning subgraph. In particular, the set E used in this section will be interpreted later as the set of 
edges incident to a given vertex. 

Let E be a finite set of elements e and Y = (Y e , eeE)e [0, oo) E . We define a probability measure 
fi on {0, 1} E as follows: the binary random variables B e with e G E are independent Bernoulli random 
variables with jj,{B e = 1) = yq-p- G [0, 1). For a given integer w, we will be mainly interested in the 
measure obtained from fj, by conditioning on ~^2 eeE B e < w and we denote it by ~p. Simple calculations 
show that for any (77(e), e e E) £ {0, 1} E such that J2 e eE ^( e ) — w > we nave 



eeE ' 

77(e) 



71 {B e = 77(e) : Ve G E) := fi [ B e = 77(e) : Ve G E 



B e <w\ = 



TleeE Ye 



eeE ) 2^SCE,\S\<w 1 



where for any set S C E, Y s is a convenient notation for Y s := Yiees^e ana - ^ = ^ an< ^ 1^1 ' s ^ ne 
number of elements in S. Note that with Y e = z for all e, we recover {1} for a simple star graph, where 
all edges have exactly one node in common with associated degree constraint w and all other nodes have 
constraints larger than one. 

For a given e G E, we introduce the following notations: 

Y e K e {Y) 



fl{B e = 1) = /i \^B e = 1 
where (we use the notation E\e for _E\{e}): 



eeE / 



TZ e (Y) := ^| S |<-iY g < 1. (4) 

Note in particular that the function Y M> 7?. e (Y) depends only on the components Yt with I ^ e. 
We also define 



P(Y) := £ M ( S e = 1 E Be < ™ ] - E r 

e££ \ eE-E / eeE 



Y e K e {Y) 

+ Y e n e (Yy 



With the interpretation given above of the star graph with all Y e = z, D(Y) is simply the mean degree 
of the central node under the probability distribution ((T|). 

The following proposition shows crucial monotonicity properties of the functions introduced above. It 
follows the argument in |26j (see Lemma 1 and 2) using results from the theory of negative dependence 
in [25]. 

Proposition 1. The mapping lZ e : [0, oo) E \ e — > (0,1] is non-increasing in each variable. The mapping 
T> : [0,oo) E — t [0, \E\) is strictly increasing in each variable. The mapping z G [0,oo) H> zlZ e {z~X.) (where 
zX denotes the vector with entries zX e ) is strictly increasing i/X > 0. 
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Proof. From Theorem 2.7 in |25) . we know that the measure [i is negatively correlated: for e 7^ /, 

Jl(B e = l,B f = l)< Jl(B e = l)Ti{B s = 1). 
This result directly gives the first point of the proposition since we have: 



dTZJY) 



E 



SC-E\{e,/}, \S\<w-2 



E 



e}, |S|<to- 



SCE\{e,/}, \S\<w-l 



dY i 



f 



E 



SC-E\{e}, \S\<w 



Y s 



(Esc-B\{e}. \S\<w 

■p(B e = Q)]l(B e = 1, B f = 1) - Ji(B e = l)ji(B e = 0, Bf = 1] 



Y s 



Y e Y f Jl(B e = 0) 2 
■p(B e = l,B f = l)- JJ{B e = lygBf = 1) 



< 0. 



Y e Y f Jl{B e = 0) 2 

For the second point, we compute the derivative of T> as follows: 

dV{Y) _ K f (Y) ^ Y e OR, 



dY f 



{l+Y f H f y ' ^ f {l + Y e n e {Y)f VY, 



f 



(Y) 



y- (P( B * = l,Bf = l)~ HBe = = 1)) 



eeE 



We need to prove that this last quantity is positive. Let \x k denote the law fi conditioned on Eee-E B e = k. 
By Theorem 2.7 in [25], this measure is still negative correlated so that for e ^ /, we have: 

JI k (B e = l,B f = 0)= JI k (B e = 1) - MBe = 1, Bj = 1) > Jl k (B e = 1)74(5/ = 0). 

In particular, we get 

JC k (B e =l\B f = 0)> JI k (B e = 1) > Ti k (B e = l\B f = 1). 



Add an extra variable -B/», so that H k {B e = l|B/» = 0) = fi \ B e = 1 E/e_e^/ = an< ^ A*fc(B e 
l\Bf* = 1) = (j, ( B e = 1 E/eB B/ = & — 1 ) ■ Hence we proved 



/i B e = 1 



2 Bf = k ] > n \ B e = 1 

feE 



/EE 



Denoting a k = ji yB e = 1;E/g_eB/ = fcj and = /i ^E/e^B/ = fcj , previous inequality shows that 
£k is non-decreasing in so that for 1 < k < w, we have S J - fc-1 ? J < < S"- J - fc ? J which gives 



/J B e = 1 



^B/^fc-l] </Z [ B e = 1 



E % ^ fc ) . 
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so that we get 



B e = 1, Y, B f > k 1 -M E B / ^ fc ^ p 2 B/ > A ) 7* B e = 1, 53 B/ < A - 1 



s e = i, 53 B f > k ] > rAj2 B f- k ] ^{Be = 1) 



feE 



JeE 



B e = l, 53B/>fc (l-7l(B e = l)) > 71 E B f >k, S e = 0)7r(B o = l) 



feE 



JeE 



E B f * k 

JeE 



B e = 1 ] > A 1 ( E B / - fc 



Be =0 



in particular, we get /x (^X1/g_e Bf > k B e = \j > fi ( X)/gB Bf >kj. Note that for fc = 0, the inequality 
is strict so that summing over A, we get 

E nBf = l\B e = 1) > E M(B/ = 1), 
feE feE 

which concludes the proof for T>. Finally, we compute: 

dzn ^ = n e{z x ) + z E x f dn ^ 



dz 



feE\e 



r(B e = i) 



E 



Jl z (B e = !,% = !)- Jl z (B e = l)jl z (Bf = 1) 



zX e Jl z {B e = 0) ' ^ zX e Jl z {B e = 0) 2 

where Jl z is the conditioned measure obtained with parameter Y = z~K. Then, we have 
dzTZ e (zX) _ 1 



dz 



zx e -i z (B e = oy 



Er(B/ = i,B e = i)-r(B f = iw(Be = i) 

feE 



It follows again from Theorem 2.7 in [25] that \i z is still negatively correlated (stability under external 
fields) so that previous calculation is still valid and this concludes the proof. □ 



In what follows, we will need to consider the extension of previous mappings in order to cover the 
case Y e = oo. We now consider Y € [0, oo] E . Let E'(Y) C E be the set of elements e such that Y e = oo. 
We now show that it is possible to extend continuously the mappings lZ e and T>. To start, it is easy to 
check that for any extension of ~p(.), we must have: 



• if |£"(Y)| > w, then YleeE'(Y) M (Be = 1) = w, so that Ve £ E'(Y), we have for any Y e > 0, 
■p(B e = 0) = 1; 

• if \E'(Y)\ < w, then for any Y e > with e £ B'(Y), we have for any (77(e), e £ B) G {0, 

71(5.= 17(e): e G B) = if ^ = ^ V ~ ~ V s ' < S > 

\e££'(Y) / ^SCE\E'(Y),\S\<W-\E'(Y)\ 1 



This simple fact allows to extend the mapping lZ e as follows: 
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Lemma 1. We define the mapping lZ e : [0, oo] B \ e — > [0, 1] defined by on [0, oo) E and on [0, oo] B \ e \[0, oo 
we let E' e (Y) = {/ 6 E\e, Yf — oo} and define lZ e by 



if\E' e (Y)\ < w, then 



SscM-EKYIulel, \s\<w-\E'nn\-i'Y 



TZ e (Y) : = ^c^t(Y)u { ,},|5|< m -| gt (Y)|-i £ (Q) x] _ 



• »/ |££(Y)| > w, then K e (Y) = 0. 
In particular, we have 



TZ e {Y) = <=> = °°)>w. (6) 

/G-B\e 

T7ie mapping lZ e is continuous. 

Proof. If \E' e (Y)\ > w, then we have for all Y e > 0, Ji(B e = 0) = 1+Ye ^ e (Y) = 1 so that ^( Y ) = °- 
If \E' e (Y)\ < w, then from © we have for all Y e 6 (0, oo), 

Y u 



EscB\E'(Y), |S|<iL.-|B'(Y)| 



C/CB\B'(Y), \U\< W -\E'(Y)\, e£U ^SCE\E'(Y), \S\<w-\E>(Y) 
J2scE\E'JY)u{e}, \S\<w-\E> c (Y) 



EscE\^(Y)u{e} ; \S\<w-\E' e (Y)\ YS + *e I^5C-E\-E^(Y)u{e}, |5|<u.-|_E^(Y)|-l ^ S 

And the lemma follows since p(-B e = 0) = 1+y ^ (y) • ^ 



It is also possible to extend the mapping T> as follows: 
Lemma 2. TTie mapping V : [0, oc] B -> [0, defined by 



W = E 1 + y^Y) 1(re<w)+WA ^ 1(re = " )) ' 

e£E e eK - ' e£E 



is continuous. 



Proof. Since lZ e {Y) € [0,1], we need only to deal with the case where there exists an element e with 
Y e = oo and lZ e (Y) = 0. In this case we have YlfeE\e ^-00 = 00 ) — w ' so * na * -^"( Y ) = {/ "= Yf = oo} 
is such that |25'(Y)| > w + 1 so that 1Zf(Y) — for all f £ E. In particular, we have 

and the lemma follows. □ 



A third operator will be important in the sequel. By Proposition [T] we can define for X £ (0, 1} E : 
Q e (X) = lim t^ e (zX) = 



Escfi\e, \S\=w-l -^- S 



SscE\e, \S\=w ^ 
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Since the mapping (2, X) G (0, 00) x (0, 1] E 1— > z7?. e (zX) is non-decreasing in 2 and non-increasing in 
X, we can extend the mapping Q e : [0, 1} E — > [0, 00] continuously by, 

Q e (X) = 00 ^ > 0) < 10. (7) 

Lemma 3. The mapping Q e : [0, 1} E — ► [0, 00] is non-increasing in X and /or any X G [0, 1] B , we have 

Proof. If Ve, either Q e (X) = 00, or X e = 0, then the left-hand side is and either there exists e such 
that Q e (X) = 00, so that by © the right-hand side also equals 0, or for all e, Q e (X) < 00 and X e = 0, 
in which case J2 e eE l(^e > 0) — 0, so that the right-hand side also equals 0. Conversely, if there exists 
e such that Q e (X) < 00 and X e > 0, then we have YlfeE ^-( X f > 0) > w + 1 so that Q/(X) < 00 for all 
/ G E. Hence we have 



X e Q e (X) 

+ X e Q e (X) 

J2\s\=w,ees -^~ S 
E 2|s|=t«X 



ST^H) 1(a(X)<TO) " i(£K*.>°>>» + i)£t 

e££ e eV ; \ee_E / eEE 

= l(^l(X e >0)>to+l) 

\e£E ) eeE 

= Wl I Y H X e > 0) > W + 1 J . 

□ 

In view of © and Q, we now introduce a fourth important operator: for any I G {0, 1} E and any 
e G E, we define 

■p e (I) = 1 // < w ) , where a sum over the empty set is equal to zero. (8) 

\feE\e J 

Clearly this operator is non-increasing in I. The local operators defined in this section will be the 
building blocks of a global operator defined on the graph G. When the graph G is a tree, there is 
a simple connection between this global operator and the distribution |T]) and we describe it in the 
next section. Before that, we introduce some notations. In the sequel, there will be more than one 
set of element E and weight w involved. Indeed a local operator will be associated to each vertex 
or (directed) edge of G. To avoid ambiguity, we denote explicitly the dependence in E, w as follows: 

-p(E,w) 5 n (E,w) ^ V{E , W ) ( Q (E,w) ^ v (E,w) _ 



3.2 Finite graphs 

We now come back to the case where G = (V,E) is a finite simple graph and w = (w v , v G V) is a 
vector of constraints. For any i G V, recall that di denotes the set of edges incident to i. We now define 
local operators with associated set of elements di, constraint Wi and some weights Y. Clearly each edge 
(ij) G E will appear in the local operators associated to di and dj with possibly different weights. In 
order to distinguish between these two cases, we will add an orientation to each edge: for the operator 
with constraint Wi , the set of elements will be the oriented edges adjacent to i and directed towards i. 
We introduce now some notations. We denote by E* the set of oriented edges. For i, j G V 2 two adjacent 



11 



nodes, we denote by i — > j the oriented edge from i to j. We still denote by di the set of oriented edges 
towards i. It will be clear from the context if we consider oriented or non-oriented edges when we write 
di. We also define di\j as the set of oriented edges towards i except j — > i. For a vector X G M. E (or 

K.^) and a subset of edges F C E (or oriented edges F C £^), we denote by X[F] the vector restricted 
to F: X[F] = (Xi, i G F). 

If the graph G is finite and acyclic, it is easy to see that the law of each marginal B[e?i] under \i z G 
is exactly described by previous local probability measures f^ dl,Wi \ with constraint: to,; set of elements: 
the oriented edges toward i, di; and for some parameters (Y-g(z), e S 9i) to be computed. We now 
describe how to compute these parameters: Y(z) = (F-^(z), ^€£)e R-^. We first define a global map 

R:Ye i-> Z G K . For a given oriented edge j — > i, we consider the local map 7^^'™^(.) as in 
previous section with set of elements di, and constraint w, so that 

9i,ti)i)/ v ^ _ SscOi\j",|S|<tu,-l 



^ Y ) 



In order to define Z = 7£g(Y), we set Z^j = 72i_^ (Y). Note that the notations are consistent with 



previous section and Z^j = TZ^^"'\y) is a function of the Y^; for £ ^ j neighbour of i. In order to 



lighten the notation, when no ambiguity is possible, we will simply denote 71^'™ by Hj. 




Proposition 2. On a finite acyclic graph G and for any z > 0, i/iere is a unique solution to the fixed 
point equation: 

Y(z) = zTZ G (Y(z)). (9) 

Moreover the marginal law of B[9i] under fi G is given by the local probability measure Jl^ l ' Wi \ with 
constraint u>i, set of elements di and parameters (Y-g(z), it G di) being the restriction to coordinates in 
di of the solution to 0). 

The fact that the recursion (O is exact on trees is a standard result in the literature on graphical 
models [23J (in this context see Lemma 3 in [22] )■ It follows directly from this proposition, that the mean 
degree of vertex v G V under [i z G defined by ^ is given (when G is a tree) by: 

V (Y(z)) V >V(*)ft*(Y(z)) 

" ( l ,) ^ l + Y^(z)n^Y(z)Y (W) 

Recall that we are interested in the limit z f oo in which case, [i z G converges to the uniform distribution 
over maximum spanning subgraphs. Iterating (j9|), we obtain Y(z) = z1Zg(z1Zg(Y (z))) . It is natural to 

define on [0,1]^, 

Q G (X) = lim t^ G (zX), 

z— >oo 

which is well defined by the monotonicity properties shown in previous section. Then, X' = <2g(X) is 
defined by: 

V - n - ^scdi\j,\s\=wi-i xS 
■^i-t-j — ^tj^i{-^) — — tFs~' 

l~<Scdi\j, \S\=Wi A 
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As will be shown latter in a more general context, we have lim^oo j" Y(z) = Y e [0, oo]^ so that by 
passing to the limit in Y(z) — zTZ G (z'R, G (y (z))) , we see that Y = Qg o 1Z g (Y). 

Recall that V e was defined in (jSJ) and we extend it to a global operator: V G ■ {0,1}^ — > {0,1}^ as 



follows: 



VI e {0,1}^, 7^(1) = 1 J2 I f <w * 

\fedi\ 3 

Proposition 3. Let I G {0, 1}^ be defined by I-g — 1(Y-^ = oo). For a finite acyclic graph G, I is the 
unique solution to the fixed point equation: I = Pg(I) an d we have 

2M(G) = hm ^ 2\(Y(z)) - £ L„l I £ > to, + 1 j + w v A ^ j . (11) 

Proof. If X = K G (Y) and I' is defined by F^ = 1(X^ > 0), then by © and ©, we see that I = V G [V) 
and I' = V G (T) so that I = V G o V G (I). Moreover, starting form the leaves and iterating V G , is is easily 
seen that V = I and that it is the unique solution to I = "Pg(I)- Then, we have 

hm5>,(Y(z)) = E[E l + yS^Y) 1( ^ <w) + W - A E **)• 
Using = ft-^(Y) and F-^ = Q_-^(X), we get 

Hm E^(YW) = E ( E l+^oSx) 1( ^ (X)<ro)+W " A E 



= E ( ^ 1 E > °) - + 1 1 + w -° A E ^ ] > 

where we used Lemma [3] to get the last equality. The proposition now follows from the observation noted 
above that 1(X^ > 0) = 1^ = 1-$. □ 

Remark 1. If we remove the assumption that the graph is acyclic, then 



• there might be no solution to the fixed point equation I = V G (T). This is the case for the complete 
graph with 3 vertices and all weights equal to 1. 

• there might exist several solutions to the fixed point equation I = V G (T). This is the case for 
the complete graph with 4 vertices and all constraints equal to 2: all solutions can be obtained by 
permuting the role of the edges in the following picture: 




Figure 2: solutions to I = V G (T). 
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If G is a finite tree, this proposition with @ shows that the size of a maximal spanning subgraph 
M(G) boils down to the computation of (fTTj) which is exactly the formula ([3]) (for the particular case 
w v = w). We will now show how to extend this computation to the case where G is possibly infinite. 

3.3 Infinite unimodular networks 

We first need to recall the general framework of [2] . We still denote by G = (V 1 E) a possibly infinite 
graph with vertex set V and undirected edge set E. We always assume that the degrees are finite, i.e. 
the graph is locally finite. A network is a graph G — (V, E) together with a complete separable metric 
space 5 called the mark space and maps from V and E* to 3. Images in 3 are called marks. Each edge is 
given two marks, one associated to each of its orientation. A rooted network (G, o) is a network with a 
distinguished vertex o of G called the root. A rooted isomorphism of rooted networks is an isomorphism 
of the underlying networks that takes the root of one to the root of the other. We do not distinguish 
between a rooted network and its isomorphism class denoted by [G, o]. Indeed, it is shown in [2] how to 
define a canonical representative of a rooted-isomorphism class. 

Let Q* denote the set of rooted isomorphism classes of rooted connected locally finite networks. Define 
a metric on by letting the distance between [Gi, ojj and [G2, 02] be 1/(1 + a) where a is the supremum 
of those r > such that there is some rooted isomorphism of the balls of graph-distance radius [rj around 
the roots of Gi such that each pair of corresponding marks has distance less than 1/r. G* is separable 
and complete in this metric [2]. 

Similarly to the space G*, we define the space £7** of isomorphism classes of locally finite connected 
networks with an ordered pair of distinguished vertices and the natural topology thereon. 

Definition 1. Let p be a probability measure on Q*. We call p unimodular if it obeys the Mass- 
Transport Principle: for Borel f : Q** — > [0, 00], we have 



Let IA denote the set of unimodular Borel probability measures on £/* . 

We now illustrate this concept with examples from our problem: a (possibly disconnected) finite graph 
G = (V, E) with a vector of constraints w = [w v , v G V) is a network with mark space R and vertices 
have marks given by the vector w (and oriented edges have mark zero) . We still denote by G the network 
and we write G x for the connected component of x in G. For a (deterministic) finite network G, we 
define XJq as the distribution on induced by the uniform measure on the vertices of G: more precisely, 
Ug (\G x ,x\) — ryT ^2 ye y 1 {{G y ,y) G [G x ,x]). It is easy to verify that XJq is unimodular for any finite 
network G since the equality (|12[) is equivalent to an interchange of (finite) summation in: 



If p is a unimodular probability measure on with marks on vertices corresponding to their degree 
constraints, it will be useful to construct other unimodular measures from p. We will always start with p 
a unimodular probability measure on with marks on vertices corresponding to their degree constraints 
(and marks on edges being zero). We will take as mark spaces: 3 = R2 so that the initial mark of vertex 
v is simply (w v , . . . ) . Note that the map TZg is still perfectly defined for any locally finite network G and 
is isomorphism invariant, in the sense that for a graph isomorphism ip, we have Ticp(G) 00 = KcfYo^ 1 ) 

where vectors in Z G M^are viewed as maps from E to R. Hence, we define 1Z\g,o] = Tic where (G 1 , of) is 
the canonical representative of [G, o]. To simplify notations, we write TZq = TI[g,o]- We define Y°(z) = 
and for k > 0, Y fc+1 (z) = zlZc(Y k (z)). It is easy to show that if we set the mark for each edge e to 




(12) 
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(Y^(z), Yi-(z), . . . , Y^.(z),0, . . . ) we should still get an unimodular probability on Q* . Indeed, write 
explicitly the dependence in the mark and define g(G, o, x, Y(z)) = f(G, o, x, zTZq(Y(z))). If we proved 
that the measure p with the marks (Y^.(z), Yh.(z), . . . , Y k f l {z), 0, . . . ) is unimodular, we get 

J xev xev 

= fY,9(G,x,o,Y h -\z))dp{[G,o]) 



f Y,f(G,x,°,Y k (z))d P ({G,o}), 
J xev 



and the claim for k follows. 



By Proposition [U we can define lim^oo t Y 2k (z) = Y (z) and lim^oo 1 Y 2fe+1 (z) = Y+(z), with 
Y+(z) = zK G (Y-(z)) and Y~(z) = zK G (Y+(z)) p-a.s. We now define for any Y G K+: 

/(G, o,x, Y) = >yo^o(Y) ^ for Q ^ ad j acentj 

J- I I X— ¥0 '**X — (, * J 

and /(G, o,x, Y) = if o and x are not adjacent. Note that we have: = Y°(z) < Y~(z) < Y + (z) < 
z = Y 1 (z), in particular all quantities are finite, hence by definition: 



Jv o (Y-(z))dp([G,o}) = Jj2f(G,o,x,Y-(z))dp([G,o}) 



Y-_+ (z)K x ^ (Y-(z)) 
S l + ^o(^ (Y-(,))^ [G '° ]) 

\- ^o^(Y+(z))y+^(z) rrr n 
IT^(«^ |G ' ° ]) 



/ £/(G,x,o,Y+(z))dp([G,o]) 
^ xev 



where we used the identity: Y x l^ (z)TZ x ^, (Y~ (z)) = 1Z ^ X (Y + (z))Y^ x (z) . Since / is continuous in Y 
and bounded by one, we can apply the dominated convergence theorem and the Mass- Transport Principle 
with marks Y 2fc+1 (z) which converge to Y + (z), to get: 

fv o (Y-(z))dp([G,o}) = J f(G,o,x,Y+(z))dp([G,o}) 

xeVG 

V (Y+(z))dp([G, o]), 

Write deg(p) for the expectation of the degree of the root with respect to p. If deg(p) < oo, then the 
above expectations are finite and since by Proposition [TJ we have p-a.s. V (Y~) < £> (Y + ), it follows 
that we have indeed equality and then by the strict monotonicity of V , we get Y 2 T_^ (z) = Y+_^ (z), for 
all x £ do and z > 0, p-a.s. Then by Lemma 2.3 [2J, this directly implies that Y~(z) = Y + (z) p-a.s. We 
denote it simply by Y(z) and we have Y(z) — zTZq(Y(z)). 

Let now p Y be a probability measure on with marks on vertices w and on edges Y(z) such that 
Y(z) = z1Zg(Y(z)) p Y -a.s. We assume that the marginal law p of p Y where marks on edges are ignored 
(and only marks on vertices w are kept) is unimodular. Since < Y(z) < z, we have by a simple 
induction Y~(z) < Y(z) < Y+(z) so that Y~(z) = Y+(z) = Y(z) p Y -a,.s. if deg(p) < oo. Hence we 
proved the following proposition: 



15 



Proposition 4. Let p be a unimodular probability measure on with marks on vertices corresponding 
to degree constraints w . If deg(p) < oo, then the fixed point equation Y (z) — z7Zq(Y(z)) admits a unique 
solution Y(z) for any z > for p-almost every (G, w). 

This proposition allows us to extend the fixed point equation of Proposition [2] to an infinite setting. 
In the case where a sequence of finite graphs G n admits a random weak limit p concentrated on trees (to 
be defined in Section H]), the solution Y(z) will allow us to define the limit of pf G . We now concentrate 
on the z — y oo limit and derive a result similar to Proposition [3] for unimodular networks. 

We first introduce a convenient definition: 

Definition 2. Let p be a unimodular probability measure on with mark space H. Let M be maps from 
ioH. We say that M is spatially invariant if the measure obtained from p by adding the marks M 
is still unimodular. With a slight abuse of notation (i.e. we still denote by p the measure with marks in 
H x £), spatial invariance is characterised by: for Borel f : x — > [0, oo], we have: 

J ]T f(G,o, x ,M)dp([G,o]) = J ]T f(G,x,o,M)dp([G,o}). (13) 

For example, we already proved that Y k (z) is spatially invariant for our generic p. 

Note that the maps Qq and Vg are well-defined for any locally finite network G and we can define 
them for the isomorphism class [G, o] thanks to its canonical representative, as we did above for TZg- We 
can now define for any I £ {0, 1}^, 

F (I) = W a l{ >Wo + l)+WoA^ 4_> . 

xGdo x£do 

Recall that for any Y(z), the 'degree' of v. V v (Y(z)) is defined by (fTU)) . We now prove the following 
result 

Proposition 5. Let p be a unimodular probability measure on with deg(p) < oo. Let Y(z) be the 
unique solution to Y(z) = z7\Lg(Y(z)) as defined in Proposition^ then Y(z) f Y as z tends to infinity, 
where Y is the minimal solution to Y — Qq o T&q(Y) and we have 



hrn j V {Y{z))dp([G, o]) = J V (Y)dp([G, o]) = inf F ( 



(I)dp([G, })>, (14) 



where the infimum is over all spatially invariant solution of 1 — Vg ° 7- > g(I)- Moreover, when p is 
concentrated on trees, we can restrict the infimum over all spatially invariant solution of I = Vg 
Vg(1) which are locally independent, i.e. such that the components of the vector I restricted to do are 
independent. 

In the case where p is concentrated on finite trees, we recover exactly Proposition [3] since in this case 
there is a unique solution to I = Vg ° ^g(I) = ^g(I)- One can check that in the example of Figure [TJ 
we have: 

f q f 14 r in 

J F (Ii)dp([G,o\) = -,j F (I 2 )dp([G,o]) = - and j F (l 3 )dp([G,o]) = 

where the last quantity is the mean degree in a maximum spanning subgraph. Note that the condition 
I = Vg ° Pg(I) is crucial for p4]) to be valid. 

In the infinite setting, there might exist different spatially invariant solutions to I = Vg°Vg(X) even in 
the case where p is concentrated on trees. In the particular case of branching processes, the computation 
of J F Q (T)dp([G, o]) for each of them is easy and we do it in the following section. In order to conclude 
and get the asymptotic size of a maximum spanning subgraph, we still need to show that we can invert 
the limits in z and n and this is done in Proposition [B] 
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Remark 2. Proposition\^is valid for any unimodular probability measure p, even if p does not concentrate 
on trees. However, we can interpret T> {Y{z)) as the degree of the root in a spanning subgraph taken at 
random according to the Gibbs measure with activity z, only when p concentrates on trees. 

Proof. For x G [0,oo], we define I(x) = l(x = oo) and I(.) acts similarly on vectors (componentwise). 
For I £ {0, 1} E , we define ooi as the vector in {0, oo} E having the same positive components as I. Note 
that we have ooi( x ) < x for any x £ [0, oo] E . 

Lemma 4. Let p be a unimodular probability measure on with deg{p) < oo. For Y spatially invariant, 
we define Y' = Q G o K G (Y). Then I(Y') = "Pq o T g (I(Y)) and 

• if Y < Y' then j V {Y)dp > J F (I(Y))dp; 

• if Y > Y' then J V Q (Y)dp < J F (I(Y))dp. 

Proof. We have 

Jv (Y)dp = / ( E ! +?°^° ( ^Y) <°°) + w ° A J2 J ( y ^°)) MP, °D 

\xddo x x x£do / 

a / ( 5 rrfeSo "^ < °°> + A 5 7<y ^ > ) * ([G ' * 

\iEOo j-^u v / xEdo / 

if Y' < Y. Let X = K G {Y) so that Y' = Q G (X) and we have thanks to the MTP applied to the first 
term: 

fv o (Y)dp > / ( E i + 'x ° Q X Q °JX) 1 (Q"->°( X ) < °°) + w ° A E J ( F ^°) J dp{[G, o]) 
= f (tool ( E l(*x->o>0)>tfl o + l) +w A E /(Y*-vo) W[G,o]). 



□ 



Let I be a spatially invariant solution of I = "Pg ° T'g(I)- We define W° = ooi and for A; > 0, 
W fc+1 = Qg o 7?-G(W fc ). By previous lemma, we have 



I(W fc+i ) = V G o P G (I(W*)) = I. 

Hence we have W° < W 1 and since both Q G and TZ G are both non-increasing the sequence W fc is 
non-decreasing and we denote by W 1 its limit. We clearly have I(W I ) > I and W 1 = Q G o 1Z G (W I ) so 
that by previous lemma: 

J V (W I )dp([G,o]) = J F (I(W I ))dp([G,o]). 

Moreover since W fc+1 > W fc , previous lemma implies 

J F (I)dp([G,o}) = J F (l(W k ))dp([G,o})> J V (W k )dp([G,o}). 
Moreover, recall that T> is increasing so that by the monotone convergence theorem, we have 

lim [ V (W k )dp([G,o]) = [ V^W^dpdCo}). 

fc->oo J J 

Hence we proved: 
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Lemma 5. Let I be a spatially invariant solution o/I = Vg ° Pc? (•"•)• Then W 1 = Qg o ^(^(W 1 ) and is 
smc/i that ICW 1 ) > I and 

J V (W I )d P ([G,o}) = J F (I(W I ))rf / ;([G,o]) < J F (I)dp([G,o}). 

Lemma 6. If Y is a spatially invariant solution to Y = Qq o 1Zg(Y), t/ien we /iaue Y = W 1 ^). 

Proof. Since W° < Y, we have W fc < Y and W 1 ^) < Y, in particular I(W I ( Y ^) < I(Y). By previous 
lemma, we also have I(W I ' Y ^ ) ) > I(Y) so that we have indeed equality and hence 

J V^W^dpdCo}) = I F (I(W I W))dp([G,o]) = J F (l(Y))dp([G,o}) = J V (Y)d P ([G,o}). 

Since T> a is increasing and W 1 ^) < Y, we finally get equality. □ 

We are now ready to prove Proposition[5] We first prove that z i— > Y -^- and z i— > Y(z) are respectively 
non-increasing and non-decreasing. We need only to prove that this result is correct for Y k for any k > 
and this is shown by a simple induction on k as follows: consider z < z'; if Y k (z) < Y k (z') then 

we have — > f — - since TZq is non-increasing by Proposition [1] if — — > — y- 1 - then we 

have Y k+1 (z) < Y fe+1 (z') since (z,X) >-> z1Zg{zX) is increasing in z and non-increasing in X still by 
Proposition [TJ 

Hence we can define lim^oo t Y(z) = Y G [0, oo]^ and ]ixa. z _ yoo 4- — X G [0, 1]-^ so that by 

passing to the limit in X^-lM = -R G {Y k {z)) and Y fc+1 (z) = zTZ G [z^^j, we have X = TZ G (Y) and 

Y = Qg(X). Again since the pointwise limit of a sequence of measurable functions is measurable, we see 
that (Y, X) is spatially invariant and Y = Qc ° 1Zg(Y). 

Let now p z be a unimodular probability measure on 5* with marks on vertices w and on edges Z 
such that Z = Q G o K G {Z) p z -a.s. The marginal p of p z where marks on edges are ignored is still 

unimodular and by definition p and p z are coupled. By definition, we have for any X G [0, 1]^ and 
z > 0, zIZg(zX-) < Qg(X) so that an easy induction on k shows that: Y 2k (z) < Z, p z -a.s. Letting first 
k and then z tend to infinity, we see that Y < Z, p z -a.s. 

Let Y° = and for k > 0, Y k+1 = Qg ° 7?<;(Y fe ) so that Y fe is non-decreasing and (by previous 
result) converges to Y, the smallest solution to Y = Qq o 7^g(Y). We have 

km J V (Y(z))dp([G,o}) = J V (Y)dp = J F (I(Y))dp([G, o]). 

Let I be a spatially invariant solution of I = T^g ° Pg (I) ■ then by Lemma [5J we have 
J J F (I(W I ))d i o([G,o]) < /" f Q (I)dp([G,o]), and since W 1 > Y, 
by Lemma |5J we get, 

J V (Y)dp([G,o}) < J V a (W l )dp([G,o]) = J FoiliW^dpdCo}). 

In the case where p is concentrated on trees, it is easily seen that Y is locally independent so that 
the last claim of the proposition follows. 

□ 
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Remark 3. With the notation introduced above, we define I k = I(Y fc ), a sequence that monotonically 
converges to I" the smallest solution to I = Vg °'Pg{^- )- We have I(Y) > I - . This implies that 
ooi- < Y so that we have Y = W 1 . Hence by Lemma\^ we get 

V (Y)dp([G,o}) = J F (I(Y))d P ([G,o}) < J F (I-)dp([G,o}). 

In particular, it is NOT necessary the case that I(Y) achieving the minimum in {1$ is the smallest 
solution o/I = Vg ° T-'g(I)- 



4 From finite graphs to unimodular trees 



In this section, we show that the analysis made in previous section on infinite unimodular networks allow 
to compute the maximum size of a spanning subgraph of a sequence of finite graphs when the number 
of vertices tend to infinity. We consider a sequence of finite networks (G n = (V n , E n , w n ))„ £ N- For 
such sequence, we will use the notion of local weak convergence introduced by [5] and [5] (see also [5]). 
Recall that we show in Section [3.31 that uniform rooting is a natural procedure to construct a unimodular 
measure on from a finite network G. This measure is denoted by Vg- For p a unimodular probability 
measure on , we say the random weak limit of G n is p if Ug„ converges weakly to p (see [Sj for details 
on weak convergence). 

Proposition 6. Let (G n — (V n , E n , w")) n6 N be a sequence of finite networks with random weak limit 
p € BIAQYVT , with \E n \ = 0(\V n \) and sup^gy^ to™ < K for a fixed K > 1. If p is concentrated on trees, 
then we have 




V {Y)dp([G, o]). 



Proof. We start be stating a simple consequence of spatial Markov property, see Lemma 4 in 26 . 

Lemma 7. For a finite simple graph G = (V, E) with a vector of degree constraints w, let v G V . If the 
induced subgraph obtained by keeping only vertices at graph distance at most 2k + 2 from v is a tree, then 
we have for any z > 0, 

V v (Y 2k (zj) <Dl=Y, (Be = 1) < V v (Y 2k+1 (zj). 

e£dv 



For a finite graph G = (V,E) and v e V, we denote by Xk(G, v) the indicator function that the 
induced subgraph obtained by keeping only vertices at graph distance at most 2k + 2 from the vertex v 
is a tree. By Lemma [7J we have 

V v (Y 2k (z)) Xk (G, v) <D*=J2»G (Be = 1) < V v {Y 2k+1 (z)) Xk (G,v) + (1 - X k(G,v))w v 

Assuming that, we have w v < K for all v, we get 

£ V v (Y 2k (z)) Xk (G,v) < -L £ Dl < _L £ V v {Y 2k +\z)) Xk {G,v) + ^ £(1 - Xk (G,v)). 

Let now G n be a sequence of finite networks with random weak limit p. For simplicity, let denote by 
p n be the probability measure Uq„ on 5* obtained by taking as root of G n a vertex in V n uniformly at 
random. We can rewrite previous inequalities as follows: 

/ V (Y 2k (z)) Xk (G n ,o)d Pn < ( D z d Pn < ( V (Y 2k+1 (z)) Xk (G n ,o)d Pn + K f (1 - Xk {G n , o)) dp n (15) 
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Note that we have V {Y k (z)) < K and by definition of the metric on the left-hand side and the 
right-hand side which depends only on the ball of (graph-distance) radius 2k + 2 around the root are 
continuous on so that we can take the limit in (|15|) . Since p is concentrated on trees, the function Xk 
is one on the support of p and we get 

J V (Y 2k (z))dp < liminf J D z dp n < limsup / D z dp n < J V (Y 2k+1 (z))dp. 

By Proposition [51 we can take the limit as k tends to infinity, to finally have for any z > 0, 

lim [ D z d Pn = [ V (Y(z))dp. (16) 



Recall that by definition, 

J \ v n\ \Vn\ 

so that we get by letting first n tend to infinity, using (|16l) and then letting z tend to infinity 
/ V Q (Y)dp = lim / V (Y(z))dp < liminf -j—M(G n ). 

J z^ca J n \V n \ 

To prove the upper bound, we start with the following observation: 



\Vn\P Gn (z) 



D z dp r , 



where the polynomial P Gn is the one appearing in the definition of Gibbs measure p Gn , see (JTJ) . Differ- 
entiating once more, we get 



z\V n \^ J D z dp n = 



d /■ zP' Gn {z) , z 2 P^(z) z 2 P' Gn (z) 2 



P Gn {z) P G Jz) P Gn {zf 

\Vn\ I J {Dlfd Pn - (J Dldp^j J > 



by Jensen's inequality. In particular, we see that z i— >• J D z dp n is non-decreasing in z > so that we 
have for any 1 < s, 

J z- 1 J D z d Pn dz < z~ x dz J D s a dp n 

Moreover, we clearly have P Gn (l) < 2> En > and P Gri (s) > s M ( G «) ; so that we get 

, v (M(G„) log s-\E n \ log 2) < J D s d Pn , 

or reorganising the terms 

\EJ log 2 



Wl M ( G n) < J D s Q dp n 



IK I logs 



Hence by letting first n tend to infinity (using the fact that \E n \ = 0(\V n \) and (jTHJ) ) and then letting z 
tend to infinity, we get 

limsup _Ljlf(G n ) < f V (Y)dp, 

n \Vn\ J 

which concludes the proof. □ 



20 



5 Computation for branching processes and asymptotics for large 
graphs 



In this section, we consider bipartite graphs G — (V, E) where V can be divided into two disjoint sets 
A and B such that every edge connects a vertex in A to one in B. We also denote by w A and w B the 
vectors of degree constraints associated to vertices in A and B respectively. We are mainly interested in 
sequences of finite networks (G n = (V n ,E n ), w A , w*)„ 6 j simply denoted by G n in the sequel, whose size 
\V n \ diverges to infinity and such that \E n \ = 0(\V n \). 

Our main assumption will be that the sequence of finite networks (G n )neN admits a random weak 
limit. Motivated by applications, we now describe a natural limit (arising when the network G n itself 
is random as described in the Introduction) which we call Bipartite Unimodular Galton Watson Trees 
(BUGWT). A BUGWT is a branching process parametrised by two distributions: a distribution tt a on 
N 2 corresponding to joint distribution of the degree of vertices in A and their associated constraints and 
a distribution ir B on N 2 corresponding to the joint distribution of the degree and mark of vertices in B. 
We will always assume that 

m A — k n k,j < 00 an d mB = ^ ^ n k,j < °°- 
k,j>i k,j>i 

In particular, the size-biased distributions associated to n A and ir B are well defined: 

~A _ ( n + ^n+l.l , ~B (" + 1 ) 7r n+l.€ 

nn - e ~ mA ' Kn - e ~ ' 

A BUGWT is then a branching process defined as follows: with probability m A l +mB , the root is of type A 
and has a random number of children and a random constraint drawn according to the distribution ir A ; 
all odd generation genitors have a distribution for the couple offspring, constraint given by n B and all 
even generation genitors this distribution is n A ; similarly with probability m A l +m ii , the root is of type B 
and has a random number of children and a random constraint drawn according to the distribution ir B ; 
all odd generation genitors have a distribution for the couple offspring, constraint given by tt a and all 
even generation genitors this distribution is tt b . The set of distributions of BUGWT is clearly a subset 
of the unimodular probability measure on and we denote it by BUQWT C hi. 

Before stating our main result in this framework, we introduce some notations: we denote by (D A , W A ) 
(resp. (N A , W A )) a random variable with distribution 7r A (resp. n A ) and by (D B , W B ) (resp. (N B ,W B )) 
a random variable with distribution ir B (resp. tt b ). For < x < 1, let D A (x) be the thinning of D A 
obtained by taking D A points and randomly and independently keeping each of them with probability x 
so that P(D A (x) =k) = Y.n>kT.^nA n k) xk ( 1 ~x) n - k . We define similarly N B (x), D B (x), N B (x). 

Theorem 3. Let (G n ) n£ N be a sequence of finite (simple) bipartite networks with random weak limit 
p G BUQWT such that D A (resp. D B ) and W A (resp. W B ) are independent with finite means. If 
\E n \ — 0(\A n \) where A n is one subset of the partition of vertices of G n , then we have 

lim J-M(G„) = M {T A {q),q = g A g B (q)}=M T A {q), 

n^oo \A n \ q£[0,l] 9€[0,1] 

where 

T A {q) = E \W A A D A (g B (q))] + ||^|e [W b 1 (D b (q) >W B + l)] , with 
g A (p) = P (N A (p) < W A ) and, g B (q) = P (N B (q) < W B ) . 

Proof. Let p A (resp. p B ) be the probability measure on obtained from p by conditioning on the event 
that the root is of type A (resp. type B): o e A (resp. o e B). The fact that 



lim -Lm(G„)= lim / V (Y(z))dp A ([G,o}), 

n^oc \A n \ z^oo J 
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follows from Proposition HJ] We will now use Proposition [5] to compute the right-hand term. 

We first adapt the expression in (fH)) to our bipartite case and express it as a function of p A . For 
x 6 [0, oo], we define I(x) = l(x = oo) and I(.) acts similarly on vectors (componentwise). By Lemma[21 
we have for any spatially invariant Y: 



£> (Y)l(o e A)dp([G, o]) 



\xedo igao / 



l(oeA)d P ([G,o}) 



We now consider a spatially invariant Y solution of Y = Qq o TZq(Y) and let X = 7£g(Y). Hence 
applying first the MTP and then using Lemma [3J we get 

£ rrgSd%o 1(y "° * 00)1(0 £ ' 4)<i ' ,,[6 ' ol) 



/ «j l ^ 1(X^ > 0) > w + 1 dp([G, o]) 
fw l(^2 Po^x(I(Y)) > Wo + l)l(o e fl)dp([G, o] 



where we used the fact 1(X X ^ > 0) = 'Po^x(I(Y)) by ©. By Proposition [5j we know that we should 
consider the minimal such Y so that dividing by J l(o e A)(ip([G, o]) = gp^qjrg] ; we obtain: 

lim / V (Y(z))dp A ([G,o}) = 



inf j /"«,„ A ^^([G^D + raj / ™ol(]T V ^ x (I)>w + l)dp B ([G,o])\, (17) 

where the infimum is over all spatially invariant solution of I = Vg°'Pg(X) which are locally independent. 
Thanks to the Markovian nature of the BUGWT, these recursions simplify into a recursive distributional 
equation that we now describe. Indeed the law of I x -*y is a Bernoulli distribution with parameter say 
q if x 6 A and p otherwise. Then the determination of p and q is done as follows: let I A ~^ B , I A ^ B be 
independent Bernoulli random variables with parameter p and I B ~^ A , lf^ A be independent Bernoulli 
random variables with parameter q; let (N A ,W A ) and (N B ,W B ) be independent random variables 
distributed according to tt a and n B respectively. Then we must have the following equality in laws: 



N J 



JV 1 



I 



ME'; 



< W J 



and, 



jB^A d_ j 



E J . 

i=l 



A^B 



< w 



Taking the expectation, we get: 

Q 
P 



\N A {p)<W A )=g A {p) 
'(N B (q)<W B )=g B (q). 



E [I A ^ B ] 
E [I B ^ A ] 

We see that F A (q) correspond exactly to the integral in ([T7)l for a solution of I = Vq ° ~Pg^) an d the 
first claim of the Proposition follows. 

We define 
F(p,q) = 



e [w A ] - : 



0<j<W A -l 



E[D A ] 



E \W B ] 



W 



E TT^d-,) 



0<j<W E 
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so that T A {q) = F(g B (q),q). 

Let 4>a (resp. 4>b) be the generating functions of D A (resp. D B ). After simple computations, we get 
(by convention, the sum over the empty set is zero): 



9 A {x) = 



h A (x) = - 



E 

<w 
(x) 



x> ct> { i +1] \l - x) 



0<j<W A -l 

dg* 
dx 



„w A -i 



*r +i) 



and 

(l-x) 



(W A - 1)! 



1(W A > 1) 



and similarly for g B and h B . 

Simple computations show that: 



OF 
dF 



E 



E 

1<j<W 

E\D A 



^l\i)g A ip) 



E[D B 



■E 



(w B -i)r B 



(1 - q)l{W B > 1) 



Hence we have 



dl 1 



^'(l)qh B (q). 



(q) = ^(l)h B (q)(q-g A og B (q)) 



dq 



Note that if P (0 < W B < D B ) = then h B {x) = for all x G [0, 1] and if P (0 < W B < D B ) > 0, then 
h B (x) > for all < x < 1. Moreover, since g A o g B (0) > and 1 > g A o g B (l), we see that (in the case 
where J- A is not constant) the local minimums of T A must satisfy q — g A o g B (q) so that 

inf [F\q),q = g A og B (q)} = inf T A {q). 
ge[o,i] ?e[o,i] 



□ 



6 Application to orientability of random hypergraphs 

We are now in position to prove Theorem [T] claimed in the Introduction. The fact that the random 
hypergraph Hn ^^h converges weakly almost surely is standard [20]. Hence we can apply our Theorem 
13] with 

D A = h, W A =£ and, D B = Poi(ch), W B = k, 
where h > I > 1 and k > 1. Note that Q(x, y) = e~ x Ylj> y fr = P(Poi(x) > y) so that we have: 

F A {q,c) = E[lABm(h,9 B ( q ))] + kQ{chq ; k + 1 \ 
g B (q) = 1-Q(chq,k), 

where we denoted explicitly the dependence in c in the expression of T A given in Theorem [3l 



23 



With the notations taken in the proof of Theorem [3J we have for h — £ > 2, 
h A {x) 



¥{Bm(h -l,x)< I) 
(h-l)...(h-£) 



h A (x) 



dh A 
dx 



(x) 



(h-l)...(h-£) 



x l ~ 2 (\ - xf- 1 - 2 (t-l + (h- 2i)x) , 



and if £ = h - 1, h A (x) = (h - l)x h - 2 and h A ' (x) = (h - l)(h - 2)x h - 3 . Similarly, we have for k > 2 
(and the convention 0! = 1), 



h B '(x) 



g B (x) = l-Q(chx,k) 
h (x) — che 
dh B 



dx 



(x) 



-chx (chx) k 1 
(k 

{chfe- chx 



(A -1)1 

(chx) k ~ 2 



(k-2)\ 



chx 



If k = 1, we have h B (x) = che~ chx and h B ' (x) = -(ch) 2 e- chx . 
We also define 

A(ar) = x - g A o 5 S (a;), 



so that we get 
A'(x) 

A"(x) 



l~h B (x)h A og B (x) 



^o 5 s (a;)(/i B (x)) «(js) 



A simple calculation shows that both and n^gp are non-increasing. Hence since <7 B is also non- 
increasing, we see that 5 is non-decreasing and A" can vanish at most once on (0, 1), so that A' admits 
at most two zeros on [0, 1] and A at most three. Moreover since A(0) = and A(l) > 0, there are 
only two cases: either A(x) > for all x G [0, 1] or there exists < u < v < 1 such that A(x) > for 
x £ [0, u] U [v, 1] and A(a;) < for x G [u, v]. In particular in this last case, u (resp. v) is a local maximum 
(resp. minimum) of T A . 

Let x*(c) be the largest solution of A(x) — 0, then we have 

min T A (q,c) = rninU, T A (x*(c), c)}, 

«6[Q,1] 

since ^^(0, c) = I and in the first case x*{c) = while in the second case x*(c) — v. A simple computation 
shows that as soon as x*(c) > 0, we have ^—(x*(c),c) < 0. Hence we need to find c such that x*(c) > 
and J rA {x*{c) 1 c) = I. Denoting £* = chx*, we can rewrite this last equation as 



cE 



(£ - Bm(h, 1 - Q(C, k))) + = kQ{C,k + 1) 



(18) 



By definition, we have £* = chP (Bm(h — 1, 1 — Q(£*, fc)) < £), so that since ^* > 0, we can eliminate the 
variable c in (TT5)) and get 



Bin(M -<?(£*,*)))" 



Q(e, fc + 1)P (Bin(/i -1,1- Q(e, fc)) < 



which is exactly the equation appearing in Theorem [T] We now show that this equation has a unique 
positive solution. By contradiction, assume that < £i < £2 are two solutions. Then applying the 
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map £ M- P (Bin(/i — 1, 1 — fc)) < ^) we can define x\ < x-i and then Cj = so that we have 

(a;i,Ci) = J 7 (x%,c%) = £ with asj = x*{ci). Since A is a non-increasing function of c, it implies 
that x*(c) is a non-decreasing function of c and ci < c 2 . But then J 7j4 (x*(c), c) = i for all c G [ci,(^], 
contradicting the relation 2^— (&*(c),c) < valid for a;*(c) > 0. 

As noted in the Introduction, H n A n\,h ^ s (A fc)-orientable if and only if M(G n ) — £\A n \. In particular, 
for c > k as defined in Theorem[TJ we see that Theorem [3] implies the result since linin^oo M y^ n ^ < £■ 
However in the case c < c* h t k , there is still a difficulty since Theorem [3] only shows that limn-^ = 
£. In words, there might exist o(n) vertices in A n with degree strictly less than £ in M(G n ). In order to 
prove that with high probability, this will not occur we need to rely on specific properties of the /i-uniform 
hypergraph H n 

This part of the proof is also done in [18] and we now describe the main steps of the end of the 
proof. First, it will be more convenient to work with a different model of random hypergraph that we 
denote H n ^ p ^'- each possible (Zi)-hyperedge is included independently with probability p. More precisely, 
given n vertices with n > h, we obtain H ntPj h by including each /i-tuple of vertices with probability p, 
independently for each of the Q) tuples. Standard arguments show that if = cn, i.e. p = chj (TZi) 
then the H n ^ p ^ model is essentially equivalent to the ifn,|cn|,ft model |20j . 

Let c < c < c* h f k and consider the bipartite graph G n obtained from H nj p t h with p — ch/{^~^\. 

Given a maximum admissible spanning subgraph S{G n ), we say that a vertex in A n (resp. in B n ) with 
degree £ (resp. k) in such a spanning subgraph is saturated. An alternating path is a path in which the 
edges alternatively are covered in S and uncovered in S. If there exists a vertex v in A n which is not 
saturated, we consider the subgraph obtained as follows: remove the edges in S adjacent to v and keep 
all alternating paths starting from v. We denote by K the resulting bipartite graph. It is easy to see that 

• the degree in K of v is at least h — I + 1 ; 

• if v G B n D K, then v is saturated in S and its degree in if is k + 1; 

• if w G A n fl K then if w is saturated in S, its degree in K is at least h — I + 1, otherwise its degree 
is at least h — I + 2. 



We now use a density argument in order to show that the size of K is linear in n w.h.p. Suppose 
that we are in a case where k > 2 (the case h — £ > 2 is similar). Let i be the number of vertices in 
B n CiK. Each such vertex is adjacent to k covered edges in S and 1 uncovered edge in S. Hence, we have 
\A n fl K\ > Hence by Lemma 4.1 of [T7], there exist 7 = 7(^1, |) such that \K C\B n \ > jn (note that 
the probability space considered in Lemma 4.1 of [TB] is slightly different from ours but asymptotically 
equivalent as shown by Lemma 3.1 of |18|). 

For c < c* h £ k , we now construct an admissible spanning subgraph on H n ,-p,h with p = chj and 
saturating all vertices in A n . Note that since p < p, there is a natural coupling between H ntPt h and 
H ntP: h' we can obtain H n>Pt h from B n ^,h by removing each hyperedge independently with probability 
p — p > 0. Start with a maximum admissible spanning subgraph M(G n ) as above. If this maximum 
spanning subgraph has all vertices in A n which are saturated, then we are done, since to obtain an 
admissible spanning subgraph on H n ^ p ^h we need only to remove the hyperedges in H n ^ p ^ which are not 
in H n ^ p ^h- If this maximum spanning subgraph has o(n) vertices in A n which are not saturated. We 
denote by gap n = £\A n \ — M(G n ) = o(n). Note that for any non saturated hyperedge, we can construct a 
subgraph K as above. If we remove an hyperedge in K , then either this edge was non-saturated and gap n 
decreased by more than one or, this edge was saturated but then at least one nonsaturaded hyperedge 
increases its degree in the spanning subgraph. So that in all cases, gap n decreased by at least one. We 
then proceed as follows: attach to each hyperedge a uniform [0, 1] random variable independently from 
each other and denoted U a , for a G A n . To obtain H UtPt h, we remove all hyperedges a such that U a < p—p. 
We can do this sequentially starting with the one with the lowest U a . Suppose that at the fc-th step (with 
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k > 1), there are still non saturated hyperedges, we consider the union of the subgraphs K described 
above. We showed that their size is at least 771. Hence with positive probability, the hyperedge removed 
will decrease the value of gap n . If it reaches zero, then we are done, otherwise, we consider the new sets 
K and proceed to the next step. Suppose that at the end of the procedure, gap n is still positive. We see 
that we decreased gap n in total by a binomial random variable with a number of trials linear in n and 
with a positive probability so that with high probability, this binomial random variable is larger than 
say en for a sufficiently small e > 0. Hence with high probability gap n reached zero by the end of the 
procedure and this concludes the proof. 
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